
###################################################################################################
### Figure 1 

canteen.1 <- c(.59, .21)
se.canteen.1 <- c(.01, .01)

par(las=1, family="Gill Sans", font.main=1, cex.main=1, cex.axis=.8, cex.lab=1)
b <- barplot(canteen.1, ylim=c(0, .6), xlim=c(.5, 3.3), border=F, 
             names.arg=c("Non Zero Sum:\nAdditional meal without pork", "Zero Sum:\n No more meals with pork"), space=.6, 
             main="Meals in Canteen Experiment",
             ylab="In favor of proposal (0-1)",
             col=c("#8da0cb", "#fc8d62"))
grid(lty=1, lwd=.5, col="lightgrey")
segments(b, canteen.1-1.96*se.canteen.1, b, canteen.1+1.96*se.canteen.1)
text(b, .05, canteen.1)


holiday.1 <- c(.40, .35, .22)
se.holiday.1 <- c(.01, .01, .01)

par(las=1, family="Gill Sans", font.main=1)
b <- barplot(holiday.1, ylim=c(0, .6), xlim=c(.5, 4.9), border=F, 
             names.arg=c("Win Win:\n Everyone gets a day off", "Non Zero Sum:\nMuslims get a day off", "Zero Sum:\nReplace a Christian holiday"), space=.6, 
             main="Public Holiday Experiment",
             ylab="In favor of proposal (0-1)",
             col=c("#66c2a5", "#8da0cb", "#fc8d62"))
grid(lty=1, lwd=.5, col="lightgrey")
segments(b, holiday.1-1.96*se.holiday.1, b, holiday.1+1.96*se.holiday.1)
text(b, .05, holiday.1)

###################################################################################################
# Figure 2

canteen.1.lo <- c(.31, .10)
se.canteen.1.lo <- c(.02, .01)
canteen.1.hi <- c(.72, .31)
se.canteen.1.hi <- c(.01, .01)

par(mgp=c(2, 1, 0), tck=-.005, cex.axis=.8, cex.lab=.8, cex.main=1)
plot(c(1,2), canteen.1.lo, type="o", ylim=c(0, 1), xlim=c(.8, 2.2), axes=F, xlab="", ylab="In favor of proposal (0-1)",
     main="Meals in Canteen Experiment")
grid(lty=1, lwd=.5, col="lightgrey")
points(c(1,2), canteen.1.hi, type="o", pch=19)
segments(c(1,2), canteen.1.lo-1.96*se.canteen.1.lo, c(1,2), canteen.1.lo+1.96*se.canteen.1.lo)
segments(c(1,2), canteen.1.hi-1.96*se.canteen.1.hi, c(1,2), canteen.1.hi+1.96*se.canteen.1.hi)
points(c(1,2), canteen.1.lo, pch=19, col="white", cex=.8)
axis(1, at=c(1, 2), labels=c("Non Zero Sum:\nAdditional meal without pork", "Zero Sum:\n No more meals with pork"))
axis(2)
box()
text(1, .15, "Low IMCP", cex=.8, pos=3)
text(1, .85, "High IMCP", cex=.8, pos=1)

holiday.1.lo <- c(.17, .17, .09)
se.holiday.1.lo <- c(.02, .02, .02)
holiday.1.hi <- c(.49, .42, .27)
se.holiday.1.hi <- c(.02, .01, .01)

par(mgp=c(2, 1, 0), tck=-.005, cex.axis=.8, cex.lab=.8, cex.main=1)
plot(c(1:3), holiday.1.lo, type="o", ylim=c(0, 1), xlim=c(.8, 3.2), axes=F, xlab="", ylab="In favor of proposal (0-1)",
     main="Public Holiday Experiment")
grid(lty=1, lwd=.5, col="lightgrey")
points(c(1:3), holiday.1.hi, type="o", pch=19)
segments(c(1:3), holiday.1.lo-1.96*se.holiday.1.lo, c(1:3), holiday.1.lo+1.96*se.holiday.1.lo)
segments(c(1:3), holiday.1.hi-1.96*se.holiday.1.hi, c(1:3), holiday.1.hi+1.96*se.holiday.1.hi)
points(c(1:3), holiday.1.lo, pch=19, col="white", cex=.8)
axis(1, at=c(1:3), labels=c("Win Win:\n Everyone gets a day off", "Non Zero Sum:\nMuslims get a day off", "Zero Sum:\nReplace a Christian holiday"))
axis(2)
box()
text(1, .02, "Low IMCP", cex=.8, pos=3)
text(1, .65, "High IMCP", cex=.8, pos=1)

###################################################################################################
# Figure 3

canteen.2 <- c(.62, .52, .40, .28)
se.canteen.2 <- c(.01, .02, .02, .02)

par(las=1, family="Gill Sans", font.main=1)
b <- barplot(canteen.2, ylim=c(0, .7), xlim=c(.5, 5.3), border=F, 
             space=c(.5, .2, .8, .2, .5), 
             main="Meals in Canteen Experiment",
             ylab="In favor of proposal (0-1)",
             col=c("#b8e186", "#f1b6da", "#b8e186", "#f1b6da"))
grid(lty=1, lwd=.5, col="lightgrey")
segments(b, canteen.2-1.96*se.canteen.2, b, canteen.2+1.96*se.canteen.2)
text(b, .05, canteen.2, cex=.8)
axis(1, at=c(1.6, 4.6), labels=c("Non Zero Sum:\nAdditional meal without meat/pork",  "Zero Sum:\n No more meals with meat/pork"), lty=0)
text(b[1], 0.68, "Climate Change", cex=.8, col="#4dac26")
text(b[2], 0.58, "Muslim Religion", cex=.8, col="#d01c8b")

holiday.2 <- c(.62, .32, .51, .24)
se.holiday.2 <- c(.01, .01, .02, .02)

par(las=1, family="Gill Sans", font.main=1)
b <- barplot(holiday.2, ylim=c(0, .7), xlim=c(.5, 5.3), border=F, 
             space=c(.5, .2, .8, .2, .5), 
             main="Public Holiday Experiment",
             ylab="In favor of proposal (0-1)",
             col=c("#b8e186", "#f1b6da", "#b8e186", "#f1b6da"))
grid(lty=1, lwd=.5, col="lightgrey")
segments(b, holiday.2-1.96*se.holiday.2, b, holiday.2+1.96*se.holiday.2)
text(b, .05, holiday.2, cex=.8)
axis(1, at=c(1.6, 4.6), labels=c("Non Zero Sum:\nEffect on economy",  "Zero Sum:\nReplace Christian holiday"), lty=0)
text(b[1], 0.68, "Womens' day", cex=.8, col="#4dac26")
text(b[2], 0.38, "Muslim holiday", cex=.8, col="#d01c8b")

###################################################################################################
# Figure 4

canteen.2.lo <- c(.52, .28, .26, .17)
se.canteen.2.lo <- c(.03, .03, .05, .04)
canteen.2.hi <- c(.67, .62, .45, .33)
se.canteen.2.hi <- c(.02, .02, .03, .03)

par(mgp=c(1.7, .7, 0), tck=-.005, cex.axis=.8, cex.lab=.8, cex.main=1)
plot(c(1:2), canteen.2.lo[c(1,3)], type="o", ylim=c(0, 1), xlim=c(.8, 4.2), axes=F, xlab="", ylab="In favor of proposal (0-1)",
     main="Meals in Canteen Experiment")
grid(lty=1, lwd=.5, col="lightgrey")
points(c(1:2), canteen.2.hi[c(1,3)], type="o", pch=19)
segments(c(1:2), canteen.2.lo[c(1,3)]-1.96*se.canteen.2.lo[c(1,3)], c(1:2), canteen.2.lo[c(1,3)]+1.96*se.canteen.2.lo[c(1,3)])
segments(c(1:2), canteen.2.hi[c(1,3)]-1.96*se.canteen.2.hi[c(1,3)], c(1:2), canteen.2.hi[c(1,3)]+1.96*se.canteen.2.hi[c(1,3)])
points(c(1:2), canteen.2.lo[c(1,3)], pch=19, col="white", cex=.8)

points(c(3:4), canteen.2.lo[c(2,4)], type="o", pch=19)
points(c(3:4), canteen.2.hi[c(2,4)], type="o", pch=19)
segments(c(3:4), canteen.2.lo[c(2,4)]-1.96*se.canteen.2.lo[c(2,4)], c(3:4), canteen.2.lo[c(2,4)]+1.96*se.canteen.2.lo[c(2,4)])
segments(c(3:4), canteen.2.hi[c(2,4)]-1.96*se.canteen.2.hi[c(2,4)], c(3:4), canteen.2.hi[c(2,4)]+1.96*se.canteen.2.hi[c(2,4)])
points(c(3:4), canteen.2.lo[c(2,4)], pch=19, col="white", cex=.8)
axis(1, at=c(1:4), labels=c("Non Zero Sum",  "Zero Sum", "Non Zero Sum",  "Zero Sum"), lty=0)
axis(2)
box()
text(1.5, .9, "Climate Change", cex=.8, pos=3)
text(3.5, .9, "Muslim Religion", cex=.8, pos=3)

text(1, .4, "Low IMCP", cex=.8, pos=1)
text(1, .8, "High IMCP", cex=.8, pos=3)

holiday.2.lo <- c(.50, .13, .40, .11)
se.holiday.2.lo <- c(.03, .02, .06, .02)
holiday.2.hi <- c(.68, .40, .54, .28)
se.holiday.2.hi <- c(.01, .02, .03, .02)

par(mgp=c(1.7, .7, 0), tck=-.005, cex.axis=.8, cex.lab=.8, cex.main=1)
plot(c(1:2), holiday.2.lo[c(1,3)], type="o", ylim=c(0, 1), xlim=c(.8, 4.2), axes=F, xlab="", ylab="In favor of proposal (0-1)",
     main="Public Holiday Experiment")
grid(lty=1, lwd=.5, col="lightgrey")
points(c(1:2), holiday.2.hi[c(1,3)], type="o", pch=19)
segments(c(1:2), holiday.2.lo[c(1,3)]-1.96*se.holiday.2.lo[c(1,3)], c(1:2), holiday.2.lo[c(1,3)]+1.96*se.holiday.2.lo[c(1,3)])
segments(c(1:2), holiday.2.hi[c(1,3)]-1.96*se.holiday.2.hi[c(1,3)], c(1:2), holiday.2.hi[c(1,3)]+1.96*se.holiday.2.hi[c(1,3)])
points(c(1:2), holiday.2.lo[c(1,3)], pch=19, col="white", cex=.8)

points(c(3:4), holiday.2.lo[c(2,4)], type="o", pch=19)
points(c(3:4), holiday.2.hi[c(2,4)], type="o", pch=19)
segments(c(3:4), holiday.2.lo[c(2,4)]-1.96*se.holiday.2.lo[c(2,4)], c(3:4), holiday.2.lo[c(2,4)]+1.96*se.holiday.2.lo[c(2,4)])
segments(c(3:4), holiday.2.hi[c(2,4)]-1.96*se.holiday.2.hi[c(2,4)], c(3:4), holiday.2.hi[c(2,4)]+1.96*se.holiday.2.hi[c(2,4)])
points(c(3:4), holiday.2.lo[c(2,4)], pch=19, col="white", cex=.8)
axis(1, at=c(1:4), labels=c("Non Zero Sum",  "Zero Sum", "Non Zero Sum",  "Zero Sum"), lty=0)
axis(2)
box()
text(1.5, .9, "Women's day", cex=.8, pos=3)
text(3.5, .9, "Muslim holiday", cex=.8, pos=3)

text(1, .4, "Low IMCP", cex=.8, pos=1)
text(1, .7, "High IMCP", cex=.8, pos=3)

###################################################################################################












